Observation of Convergent Oscillations of the Flux Line Lattice as a Result of 
Magnetic Flux Jumping in Hard Superconductors 
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We have monitored new peculiarities of the dynamics of catastrophic avalanches of the magnetic 
flux in superconducting Nb, Nb-Ti, and YBaCuO samples: i) convergent oscillations of the magnetic 
flux; ii) a threshold for entering the huge flux avalanches in the shielding experiments; iii) a threshold 
for the exit of a residual flux in the trapping experiments. The observed phenomena are interpreted 
in terms of a theoretical model which takes into account the inertial properties of the vortex matter. 

PACS numbers: 74.60.Ec; 74.60.Ge 
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Recent studies of the marginally stable Bean critical 
state determined by a balance between intervortex re- 
pulsion and pinning by defects in type II superconductors 
have brought out a number of astonishing phenomena as- 
sociated with spatiotemporal dynamic properties, vortex 
avalanches of many scales, self-organized criticality and 
memory effects in vortex matter (for example, |IJ). 

Our experiments identify flux avalanches as nonlocal 
processes propagating rapidly over the distance compa- 
rable with the size of the sample. Earlier thermal effects 
in massive rods of niobium caused by giant flux jumps 
were observed in [B. Some observations of the behav- 
ior of the flux jumps and its structure were presented in 
H . These data show that there exists the exact correla- 
tion between the jumps in the resistivity and the collapse 
of magnetizaton. One of the recent investigations (f| of 
magnetic field screening by the superconductors in the 
vortex state revealed "negative vortices" whose penetra- 
tion leads to the expulsion of a magnetic field. 

In this paper we report about the observation of two 
distinct magnetic phenomena in the vortex state of su- 
perconductors: i) a threshold for entering the huge flux 
avalanches in the shielding experiments (in other words, 
the increase in the screening properties of a superconduc- 
tor before the flux jump). We have also found a thresh- 
old for the exit of a residual flux in the trapping exper- 
iments; ii) convergent oscillations of the magnetic flux 
arising from flux jumping. 

The latter result is of a particular interest since collec- 
tive modes in a fluxoid system are very difficult to ob- 
serve IEJ. In addition, due to a large vortex viscosity, the 
displacement waves in a vortex lattice do not propagate. 



The observed oscillations in the vortex system can be 
interpreted as a result of the existence of a definite value 
of the effective vortex mass, i.e. this phenomenon can 
be considered as displaying the inertial property of the 
vortex matter. 

We have studied the structure of the flux jumps by 
means of a miniature Hall probe sensor. It was placed 
in the center of the sample (Fig. 1) and measured the 
surface induction Bq. We have examined the Hall sensor 
voltage directly by the transient recorder (model TCC- 
1000, Riken Denshi Co., LTD). In our experiment a flux 
avalanches are detected in real time (t rea i > 10~ 6 s) with- 
out distortion. In the shielding experiments, the zero 
field cool mode (ZFC) was used. The trap field mode was 
implemented by increasing the external magnetic field up 
to the values exceeding the critical magnetic field H C 2- 
The data presented in this study were obtained from Nb- 
Ti(at% 50) samples, polycrystallinc Nb plates, and melt- 
textured YBaCuO slabs. The sizes (Lxlx2Rmm 3 , 
see Fig. 1) of the samples were as follows: for Nb-Ti, 
15 x 2.8 x 3.1mm 3 ; for Nb, 11 x 3.5 x 5mm 3 ; for YBaCuO, 
6 x 2 x 3.5 mm 3 . 

Fig. l(a-c) shows the temperature and temporal de- 
pendence of the surface magnetic induction Bo for the 
first and the second flux jumps in a shielding mode for 
the Nb-Ti sample. The structure of the jump essentially 
differs from a simple step. Some of the principal differ- 
ences are clearly seen in Fig. 2(a). 

The development of the thermomagnetic instability 
can be divided into three stages (see Fig. 1(b) and 
Fig. 2(a)). The flux jump is preceded by a certain phe- 
nomenon under which the field on a superconductor sur- 
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FIG. 1: Above - geometry of the experiment. Main frames: 
the temporal evolution of the surface magnetic induction Bo 
for the first flux jump (a, b) and for the second flux jump (c) . 
Inserts: Hysteresis loops, B(H) 



face decreases (see Fig. 2(a) for Nb-Ti, Fig. 3(a) for Nb, 
and Fig. 4 for YBaCuO). The value of the negative peak 
in the experiment with Nb is 16% of the total value of the 
flux jump. This is the first stage of the process. The sec- 
ond stage is in fact the avalanche-like penetration of the 
flux. Finally, at the last stage the relaxation of thermal 
and conductive properties of the superconductor occurs. 
In the most cases, this is a nonmonotonic process. Let 
us analyze it for the Nb-Ti sample. 

As is shown in Fig. 2(a), after the avalanche had pen- 
etrated the sample, the induction Bo(t) in the center of 
the sample was higher than that of the external mag- 
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FIG. 2: (a) The fine structure of the first jump for the Nb-Ti 
sample; (b) The results of simulation of the oscillation process 
in the vortex system; insert: distribution B vs distance x for 
the first and the second jumps 



netic field by AB comp . Recently it was found by Nowak 
Q , that the large avalanches are most likely system span- 
ning; as a result after avalanche a significant fraction of 
the Nb film has magnetic induction large than the applied 
field. Apparently, Coffey JjJ observed such a compression 
of the vortex lattice by measuring the distribution of the 
induction inside the sample after the flux jump. Next, in 
the background of a relaxation curve the damped oscil- 
lations with a characteristic frequency v ~ 2.5 kHz were 
observed. 
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FIG. 3: Structure of the avalanche jumps for (a) shielding 
and (b) trapping regimes for Nb. Inserts show Bo{H). 



The effect of the compression of the magnetic flux 
caused by avalanche passing is particularly pronounced 
in the YBaCuO samples as well. In order to check this 
effect we measured the difference B — fj, H using two 
Hall probes where H is the external magnetic field, (iq is 
the permeability of vacuum. Fig. 4 shows that this differ- 
ence changes its sign after the avalanche has passed. The 
difference AB comp is about the value of the penetration 
field H p . 

Now, let us consider peculiarities of the structure of 
the thermomagnetic avalanche for the case of magnetic 
flux trapping. For Nb, (Fig. 3(b)) a positive peak Bo(t) 
is observed which precedes the development of the insta- 
bility. It is of the same origin as the negative peak Bq (£) 
in the case of the shielding mode (Fig. 3(a)). 

The above-described peculiarities of the flux jumps are 
not sensitive to the change in the rate dH/dt of the exter- 
nal magnetic field sweeping in the range of 0.05 — 2 T/s. 

Now we would like to scrutinize the observed phe- 
nomena. In our opinion the most striking effect is the 
compression of the magnetic flux after the avalanche has 
passed through the sample. It is very important to em- 
phasize that the magnetic induction in the middle of the 
sample can exceed the value of the external magnetic 
field. In this case, some of the vortices move in the di- 
rection from the area with low vortex density towards 
the area with the vortex surplus, against the total repul- 
sion force. This means that these vortices move in spite 
of all the forces acting in the direction opposite to this 
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FIG. 4: The temporal dependence of the difference of mag- 
netic induction in the middle of the sample Bo and the exter- 
nal field {1qH for the YBaCuO plate. 

motion. Such a motion is possible only due to the exis- 
tence of some effective vortex mass. There are plenty of 
papers dealing with the theoretical attempts to estimate 
the mass of the Abrikosov vortex @, |J. Experimental 
data on microwave field attenuation contain information 
on the inertial properties of a vortex Q . 

To explain the observed oscillatory phenomena, we 
suggest the following mechanism which takes into ac- 
count the inertial properties of the vortex system. Prior 
to the jump, the mixed state of superconductors is char- 
acterized by nonuniformly distributed magnetic induc- 
tion localized near the surface. As a result of the 
avalanche, the flux rushes from either sides of the sample 
towards the center. Two fronts of the penetrating flux 
collide in the center of the sample and, owing to the ex- 
isting vortex mass, give rise to the local surplus density of 
the magnetic flux that exceeds the value of the external 
magnetic field. The repulsion force in the vortex struc- 
ture at the center of the sample that have resulted from 
its compaction, initiates the wave of the vortex density 
of the inverse direction of propagation. Upon reaching 
the surface, this wave is reflected from it. This results 
in the oscillations in the vortex system. The limitation 
of the number of oscillations observed is caused by the 
existence of damping. One succeeds in observing the os- 
cillation of the vortex density only owing to a strong com- 
pression of the vortex structure as a result of the giant 
avalanche-flux. But for the second jump (Fig. 1(c)), a 
smooth descending slope in the curve Bq (t) is noted and 
the oscillatory mode is not observed. 

The observed process of damping oscillations can be 
described in terms of the following theoretical model. We 
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use the continuity equation for the dimensionless vari- 
ables, 

b+{bv)' = (1) 

where b = B/H p — Bc/2irJ c (0)d is the normalized mag- 
netic induction, d is the sample thickness, J c (0) is the 
critical current density at the initial temperature (before 
the avalanche has passed), the dot denotes the derivative 
with respect to the time r = t/to = t&oH p /4:ird 2 r], $o is 
the flux quantum, r\ is the viscosity coefficient, the sign 
prime denotes the derivative with respect to the coordi- 
nate £ = x/d, v — V/Vo — Vto/d is the hydrodynamic 
vortex velocity. The second equation of the model de- 
scribes the motion of the vortex system under the action 
of the Lorentz, pinning, and viscosity forces, 

fj,v'v + fiv = -b' - 2sign(v)f(0) - v. (2) 

Here /i = m/mo = m$>oH p /4ird 2 ri 2 is the normalized 
vortex mass, 9 = AirCAT/H 2 , C is the specific heat, 
AT = T — To is the deviation of the temperature from 
its initial value To before the flux jump, the function f(9) 
describes the temperature dependence of the critical cur- 
rent, J c (6) = J c (O)/(0), f(0) = 1 at 9 = and f(9) = 
at 9 = T c — T) (T c is the critical temperature). The 
temperature change occurs due to the Joule losses. Ne- 
glecting the process of thermal conductivity, we can use 
the following equation of the energy balance: 

9 = bv 2 + 2\v\bf(9). (3) 

This set of equations should be solved within the spa- 
tial interval —1/2 < £ < 1/2 together with the boundary 
and initial conditions, 6(1/2, t) = b(— 1/2, r) = H/H p = 
a, b((,0) = a(l + 2(/?A*)(|£| - 1/2)), = 
0, v(£,0) = 0. Here (3 = Jq/J c , Jq is the critical cur- 
rent density at the initial time moment. We assume that 
this value exceeds slightly the value J c (0), i.e. /3 > 1. 
This supposition is introduced to initiate the onset of an 
avalanche. 

The results of the simulation for a = 1, (3 = 1.01, 
f(9) = 1 - k6, k = H 2 /At:C(T c - T ) = 12 are presented 
in Fig. 2(b). It is seen that the curves in this figure de- 
scribe qualitatively well the behavior of the experimental 
ones (see Fig. 2(a)). 

The suggested mechanism of the oscillations which 
takes into account the inertial properties of the vortices 
allows one to interpret the surprising behavior of the vor- 
tex system just before the avalanche passage when the 
magnetic flux near the middle of the sample unexpect- 
edly moves in the opposite-to-avalanche direction (see 



Figs. 3(a), (b)). The mechanism of this effect may be as 
follows. Due to the vortex mass, the frozen profile of the 
magnetic induction after some avalanche has passed can 
contain a smooth peak in the middle of the sample (see, 
for example, Q). This local maximum remains on hand 
while the external field increases and exists in the initial 
profile for the subsequent avalanche. The new avalanche 
can destroy this peak due to the heating effect before 
the new vortices arrive there. This tends to decrease the 
vortex density in the middle of the sample. 

Finally, we remark that very interesting oscillating 
phenomena originated by the thermomagnetic avalanches 
are observed in the vortex matter for both low-T c and 
high-T c hard superconductors. These phenomena can 
be qualitatively interpreted in terms of the theoretical 
model which takes into account the inertial properties of 
the vortex system. However, new experimental and the- 
oretical in-depth studies are necessary to elucidate the 
situation. 
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